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Abstract 

We wish to tile a rectangle or a torus with only vertical and horizontal bars of a given length, such 
that the number of bars in every column and row equals given numbers. We present results for particular 
instances and for a more general problem, while leaving open the initial problem. 
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1 Introduction 



In general terms, tomography is the area of reconstructing objects from lower dimensional projections. We 
consider the problem of reconstructing a rectangular grid from projections on the columns and on the rows. 
Think of the grid as a layer in a crystal, and in order to measure it, we send beams through the crystal from 
two orthogonal directions. Measurements will give us quantitative information about columns or rows of the 
grid (H). Consider the problem, where each cell of the grid (think of it as an atom) is matched with at most 
one of its immediate neighbors (think of it as a chemical connection). Physicists call them monomer- dimer 
systems. Many research has been done about counting the number of configurations of such a system (|^), 
or about almost uniform randomly sampling configurations (Q). 

We are interested in the particular problem, were each cell is matched to exactly one neighboring cell. 
These objects correspond to domino tilings of the grid. A measurement will reveal the number of vertical 
dominoes in each column and the number of horizontal dominoes in each row. Given these numbers we wish 



to reconstruct the grid, or any grid which satisfies the projection constraints. As a 
natural generalization of this combinatorial problem we are interested in the tiling 
of the grid with horizontal bars of length h and vertical bars of length v, for some 
integers h,v. We call it the Tiling with Bars Reconstruction problem. Given 
a pair of column and row vectors (to, n) (tomographic constraints) and integers h,v 
we want to construct a tiling with bars, such that to counts the number of vertical 
bars in the columns and n counts the number of horizontal bars in the rows. This 
problem has two variants, whether we tile a rectangle or a torus. 
The problem is left open by this paper, but we were able to find solutions for sub- 
problems and for a more general problem. We summarize our results in this table: 
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The quadratic algorithm has been found independently by Picouleau (|3|) for a more general condition 
(see end of section pO| ). 



2 Definitions 

Let a,b > 1. The rectangle R a xb is the product {0, ■ • ■ , a — 1} x {0, • • • , b — 1} and the torus T aX b is the 
product Z a x Zb. Columns are numbered from left to right and rows from top to bottom (see Figure |^). 
Each element of a rectangle or a torus is called a cell. 
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Fi gure 2: A rectangle R a xb and a torus T a xb- 



Let /i, t> > 1. The horizontal bar of length h is the rectangle Rhxi and the vertical bar of length v is the 
rectangle R\xv If the length is 2 we call the bar a domino. 

A rectangle i? QX 6 (respectively a torus T ax &) is said to be tillable with the vertical and horizontal bars 
(of lengths v and h respectively) if it can be partitioned into those bars. The projections of such a tiling is 
the pair of vectors (m, n) = (mi • • • m a , ri\ ■ ■ ■ rib) 6 N° x N fc such that for every column i, rrii is the number 
of vertical bars in it, and for every row j, rij is the number of horizontal bars in it. 

We define the following reconstruction problems. 

Tiling a Rectangle (respectively Torus) with Bars under Tomographic Constraints 
input (m, n) 6 N° x N' and h,v>l. 

output a tiling of the R a xb (respectively T a xb) with projections (m, n). 



3 Uniform constraints 

In this section we characterize valid instances for the special case when the constraints vectors are uniform, 
that is Vi : rrii = to, Vj : rij = n for some integers to, n. Both the torus and the rectangle case are studied. 



3.1 The torus case 

The number of cells covered with vertical bars is amv and the number of cells tiled with horizontal bars 
bnh. Clearly these numbers must add up to the total number of cells ab. In this section we show that this 
condition is sufficient for a torus tiling to exist. 

Lemma 1 // a, b, h^v^m^n > 1 are such that ab = amv + bnh then there exist p, q,a' ,6' > 1 satisfying 

• (P + ?) = gcd(a, b) and a = (p + q)a' and b = (p + q)b' . 

• nh = pa' and mv = qb' . 
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Proof. If we denote c = gcd(et, b), a = ca' , and b = cb' , then the equality ab — amv + bnh can be rewritten as 
ca'b' = a'mv + b'nh. It follows that (ca' — nh)b' = a'mv. From Gauss theorem, a'\(ca' — nh), and therefore 
a'\nh. In other words, there exists p such that pa' = nh. By symmetry, there exists q such that qb' = mv. 
Finally, notice that 

cicb'nh + ca! mv) 



Theorem 1 Let a, b, h, v,m,n > 1. Let (m, ri) = (m ■ • -m, n- ■ -n) G N a x N b . The torus T a xb is (to, n)- 
tillable with the bars Rhxi and Rixv if and only if ab = amv + bnh. 

Proof. If we assume that the torus T aX b admits an (to, n)-tiling with the bars Rhxi and Rixv then, by 
simple considering the area covered by the tiling, it is direct to notice that ab = amv + bnh. 
Conversely, if ab = amv + bnh then, by Lemma [j], there exist p, q, a' , b' > 1 such that: 

• (p + 0) = gcd(a, b) and a — (p + q)a' and b = (p + q)b' . 

• nh = pa' and mv = qb' . 

As it appears in Figure ||, the torus T aX b can be partitioned into (p + q) 2 rectangles defined for each 
i, j G {0, • • • ,p + q — 1} as follows: 

0;.j = the rectangle R a >xb' whose upper left corner is the cell (a'i,b'j). 

Let us define for each i G {0, ■ • • , p + q — 1}, the following rectangular regions of T ax i,: 

• T~U — Ufct^ 1 which is simply the rectangle R n hxb' whose upper left corner is the cell (a'i, b'i), 

• Vj = Ui2+i ®i,fci which is simply the rectangle R a >xmv whose upper left corner is the cell (a'i, b'(i + l)). 

It is easy to notice that every 0j ;J - belongs to exactly one of the rectangles {"Hi, V;}o<i<p+ 9 and that therefore 
the latter is a partition of the torus T ax i,. 

In order to conclude, notice that each TLi is tillable by using only horizontal bars i?; iX i with each row 
having n bars. In the same way, each Vj is tillable by using only vertical bars Ri xv with each column having 
to bars. □ 
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Figure 3: A (2 • • • 2, 3 • • • 3)-tiling of Ti 5x i by R 2xl and i? iX 3- 



Corollary 1 If a torus T aX f, admits a tiling with uniform tomographic constraints then gcd(a, b) > 1. 
Proof. From Theorem [l] together with Lemma [l], gcd(a, b) = p + q with p, q > 1. □ 
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3.2 The rectangle case 

Theorem 2 Let a,b,h,v,m,n > 1. Let (m, n) = (m • • • m,n • • • n) € N a x N b . 77ie rectangle R a xb is 
(to, n) -tillable with the bars Rhxi &nd R\ xv if ond only if ab — amv + bnh, h\a, and v\b. 

Proof. Let us suppose that the rectangle R a xb admits an (to, n)-tiling with the bars Rhxi and Ri xv - By 
simple area considerations, it holds that ab — amv + bnh. The fact h\a follows from this observation. Since in 
every column we have mv cells tiled by vertical bars, the remaining k = b — mv cells are tiled with horizontal 
bars. Therefore k horizontal bars are between column 1 and column h, another k bars between columns 
h + 1 and 2h, and so on. In the same way we conclude v\b. 

For the converse, let ab — amv + bnh, h\a, and v\b. We reduce this case to a 01-Matrix Reconstruc- 
tion Problem. Let p, q be such that a — ph and b — qv. Now R aX b may be partitioned into pq rectangles 
Rhxv and each of these rectangles may be tiled by using one type of bars (vertical or horizontal). We define 
a p x q 01-matrix, where each entry corresponds to a Rhxv rectangle, and contains "1" if the later is tiled 
with vertical bars, and "0" otherwise. The problem is reduced to the following: given p,q,m,n > 1 such 
that pq — pm + qn, construct a 0-1 matrix of size p x q in such a way that each column has to l's and each 
row has n 0's. The solution is trivial. In fact, it suffices to consider the l's as vertical bars of unitary length 
(simple squares), the 0's as horizontal bars of unitary length (simple squares), and to apply Theorem |l| (for 
unitary length bars the torus is equivalent to the rectangle). □ 

4 Horizontal bars of unit length 

In this section, we give a polynomial time algorithm for reconstructing a rectangle tiling with horizontal bars 
of unit length. We assume in this section that h = 1 . For technical reasons we will even give a more general 
algorithm for reconstructing tilings of histograms. 

Definition 1 A histogram H of a rectangle R a xb is a subset of R aX b such that if cell is an element 

of H with j < b — 1, then (i,j + 1) is also an element of H. The top of column i is the cell (i,j) G H 
with minimal j ( remember rows are numbered from top to bottom). The number of cells of the row j of H is 
denoted by Cj . The height of column i of H is the number of cells in it. 

We will give an algorithm which, given a vector (too, mi, . . . , to q _i, no, m, . . . , rtb-i) of integer coordi- 
nates, an integer v and a histogram H of rectangle R a xb, constructs a tiling of H with the bars i?ixi and 
Rixv satisfying the tomographic constraints (to, n) or answers "No" if there is no such tiling. 

4.1 Algorithm 

This algorithm is based on a very simple idea: A solution is constructed iteratively row by row, where 
vertical tiles are placed in columns of largest remaining constraint. See http://www.lri.fr/~durr/VertOnly/- 
vertOnly.html for an implementation. 

input: to e N a ,n 6 N b ,v > 0, histogram H C R a xb- 
promise: £^=0 vm l + Y^Zl n 7 - = 
For j* from to b — 1 do 

If Cj- < nj* answer "No" and stop. 
While Cj- > nj, do 

If j* + v — 1 > b — 1 answer "No" and stop. 
Let i be a column with maximal to^ which satisfies (i,j*) G H. 
Place a vertical bar between (i,j*) and (i, j* +v — 1). 
For k from j* to j* + v — 1 do 
Remove cell (i, k) from H. 
Update Cfe := Cfc — 1. 
Update mi := m,i — 1. 
While 7J,» > do 
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Let i be any column which satisfies S H. 

Place a horizontal bar (cell actually) on 
Remove this cell from H . 
Update Jij* := rij* — 1. 

4.2 Analysis 

The correctness of the algorithm is a consequence of the lemma below: 

Lemma 2 Let T be a tiling of an histogram H satisfying the constraints given by vector (mo, ■ ■ ■ , m a -i, no, ... , rib—i). 
Let Saigo be the set of columns of maximal height where vertical bars are placed in the first step of the algo- 
rithm. 

Assume that there exists a solution. For each tiling T which solves our problem, let St be the set of 
columns with maximal height in H whose top is covered by a vertical bar. Let T* be such a solution such 
that Saigo H St* is maximal. Then S a i go = St* 

Proof by contradiction. Since |<S Zgo| = \St* | there is a column i\ of S a i go which is not an element of St*, 
and a column 12 of St* which is not an element of S a i go - Fix such a pair of columns i\, Notice that 
m l2 < m h . 

Let V\ (respectively V2) denote the set of vertical bars of T* included in column i\ (respectively column 
i-i). We construct a bipartite undirected graph G whose set of vertices is Vi U Vz. Two tiles are joined by an 
edge if they cross a same row. Notice that, in G, each bar has at most two neighbors, and G has no cycle. 
Hence, G is formed from disjoint chains (each bar with no neighbor is considered as a chain of null length). 
For every chain C let Iq be the set of rows traversing the bars in C. Then clearly Lc is an interval, and 
different chains have disjoint row sets. 

From such a chain C, one can construct a tiling Tq transforming T* , by ij i 2 q i 2 

an exchange on chains, which is the exchange of lc X {i\ } with Iq x {12}- 
Notice that this operation preserves the tomographic constraints on the 
rows, while preserving those on the columns if and only if C has an even 
number of vertices. 

Let Co be the chain with lowest row indices. 

If Co has an even number of vertices, then the tiling Tc contradicts 

the maximality of the intersection S a i go H S T *, which, consequently, Figure 4: The exchange of I c x {h} with 
achieves the proof. Ic x ^ %2 ^' 

If Co has an odd number of vertices (i.e. both endpoints of Co are bars in column 12), then by the inequality 
TOj 2 < TOjj there exists another chain C\ with an odd number of vertices, whose extremities are bars in 
column i\. Let T" be the tiling obtained from T* by exchanging Co and C\. T' satisfies the same vertical 
and horizontal constraints as T* , and, consequently contradicts the maximality of the intersection S a i go (^ST* ■ 
This last fact achieves the proof. 

□ 

Lemma 3 The algorithm presented in this section gives a tiling satisfying the constraints, if such a tiling 
exists. It's running time is 0(a log a + ab). 

Proof. We prove its correctness by induction on the number of cells of the histogram H given as input. If 
H is empty, the result is obvious. 

Now assume that the theorem holds for each histogram which has less cells than H . If H admits a 
tiling with constraints, then, by the previous lemma, there exists such a tiling using tiles placed in the first 
execution of the loop. 

After the first execution of the loop (and updating), we have to prove the theorem for an histogram which 
has less cells than H , which is true by induction hypothesis. 

Now we turn to the proof of the time complexity. The algorithm will maintain an ordering on n. 

• the initialization costs 0(a log a) time units (because of the ordering of the columns), 
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Figure 5: Trace of the reconstructing algorithm on Example [j]. 



• each passage through the loop costs 0(a) time units, since the update of the order can easily be done 
in 0(a) time units, 



• there are b passages through the loop. 
This proves that the execution of this algorithm costs 0(a log a + ab) time units. 



□ 



Example 1 For the rectangle Rioxs, for vertical bars of size 2 and the •projections 

m = (2, 3, 3, 1, 3, 3, 4, 3, 2, 2) E N 10 and n = (4, 2, 3, 3, 4, 3, 2, 6) £ N 8 . 
Figure shows the trace of the algorithm. Numbers in the cells indicate the remaining column constraints. 

4.3 Application 

The previous algorithm can be used to reconstruct a tiling, when some particular promise on m is given. 
This promise is fulfilled in particular if m is uniform, or monotone (mo < .. . < m a _i), as shown in (|5|). 

Theorem 3 Let a,b,h,v > 1. Let (m,n) = (mo • • • m a _i,no ■ • • nb-i) E N a x N b with mi — rrij for all 
i,j E {0, . . . , a — 1} satisfying [i/h\ = \_j/h\ . There is an algorithm in 0(aloga + ab) that decides whether 
the rectangle R a xb is (m,n) -tillable with the bars Rhxi and Ri xv and if yes outputs a valid tiling. 

Proof. By the same argument of the first part of Theorem^ it can be concluded that the tiling of R aX b may 
be partitioned into § tilings of rectangles of type Rhxb- It suffices now to divide every horizontal measure by 
h (i.e, to change the horizontal scale) in order to reduce the original problem to a new one in which a' = % > 
b' =b, h! = 4 = 1, v' = v, (m)' = (m ,m h , . . . , m( /_i)/j) E N a ', (n)' = nEN b (see Figure |). We can apply 
now Lemma a. □ 



5 Tiling a sub-grid 

In the previous section we showed that some instances of the TILING with Bars Reconstruction problem 
have a polynomial solution. In this section we show that a more general problem is NP-hard. In the Sub- 
grid Domino Tiling Reconstruction Problem From Projections we are given a only sub-grid 
S C R aX b to tile. We show that this problem is NP-hard by a reduction from the following problem. 

The 3-Color Consistency Problem 

We fix a set of colors A = {colorless, red, blue, green}. 



G 




Figure 6: Reducing the problem by changing its horizontal scale, 
input m c e N a and rf £ N b for every c E A. 

decide if there is a matrix T € A oxb with projections (m c ,n c ) c eA that is for all colors c we have 

< = \{j ■ Tij = c}\ and n) = \{i : T {j = c}\. 

It has been shown in ([!]) that this problem is NP-hard in the strong sense, while the 1-Color Consistency 
Problem is solvable in linear time @. The 2-Color Consistency Problem is still open. 

5.1 Gadget 

A sub-grid is a cycle if every cell has exactly two (horizontal or vertical) adjacent neighbors, and if every 
pair of cells is connected by transitivity. We start by giving some facts about cycles. 

The cycle length is always even. This can be easily seen by coloring the cells checkerboard wise black 
and white. Then adjacent cells have different colors. The claim follows from the fact that the cycle is closed. 

There are exactly two domino tilings of a cycle. Fix any numbering of the cells such that every cell 
i has neighbors i — 1 and i + 1 modulo the length of the cycle. Then clearly one tiling covers all pair of cells 
(2i, 2i + 1) with a domino, while the other one covers (2i, 2i — 1) for all i. 

We specify now a sub-grid S consisting of two cycles intersecting at a corner and an addition cell. This 
additional cell must, in a domino tiling, be matched to a cell of one of the cycles, therefore "forcing" it to 
admit a unique tiling, while the other one admits the usual two tilings. As a result we will have exactly four 
tilings of S. We define S to be the subgrid shown with its tilings in figure 0. We refer to these tilings as 
Tcoioriess,T re d,Tbi ue ,Tg reen respectively. Let (s°,i c ) be their projection vectors for every color c. Note that 
by the symmetry of S we have {s c } c <eA = {£ c }ceA- 



"colorless" "red" "blue" "green" 

2003101 1002101 2002001 2003002 




Figure 7: All four tilings of 5. 



Lemma 4 The vectors {s c } c£ a o,re linear independent. 
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Proof. Let u € N 7 be an arbitrary linear composition of the column vectors. We have to show that the 
coefficients in u = ^2 C a c s c are uniquely defined. We have 
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This system of equations has a unique solution which concludes the proof: ot C oiorie 

2 ? Gtgreen !• 



-2, C^red — 3, Qblue 



□ 



5.2 The proof of NP-hardness 
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Figure 8: Idea of the reduction 



Theorem 4 The Domino Sub-grid Tiling Reconstruction Problem From Projections is NP- 
hard in the strong sense. 



Proof. Let / = {m c ,n c ) an arbitrary instance of the 3-Color Consistency problem for an a x b matrix. 
We construct an instance of the former problem, such that there is a bijection between the respective sets 
of solutions. This proves then the theorem, (see figure 0) 
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Figure 9: Separations for h = 2, v = 2. 



We define the instance I' as (m, n) £ N 7a x N with 



m 



^2 m t s C J and " = (X) ^ 



i=l \ c / 

where ® denotes the concatenation of vectors and 

a-l 6-1 

s ' = U U^+( 71 '^'))- 

In a tiling of 5" every 7x7 block contains one of the four tilings of S. Therefore there is a natural 
bijection / between the set of domino tilings of S and the set of matrices A axb . It follows from lemma ^ 
that T is a solution to the instance I' if and only if /(T) is a solution of /. Moreover for the unary encoding 
the size of f(T) is linear in the size of T . □ 



6 Concluding remarks 

We will conclude with a observation for the general problem. Let (m, n, h, v) be an instance for the recon- 
struction problem for tilings of an a x 6-rectangle. 

We say that a particular realization is separated between column i and column i + 1, if there is no 
horizontal bar traversing the border in between. We claim that if this is the case for one realization, it holds 
for all other realizations as well: Let Ci be the number of horizontal bars beginning in column i and ending 
in column i + h — 1. Let us also denote Ci — for all i < 0. Then clearly the following induction holds: 

Ci = b - (vrrii + Ci-h+i + ■ ■ ■ + c*_i). 

Therefore there is a separation between column i and column i + 1 if and only if Ci-h+2 + • • ■ + Q = 0, which 
is a realization independent condition. 

In the same manner we define the separation between lines. These separations, which can be computed 
in linear time, partition the grid into separated rectangles which are surrounded either by a separation or by 
the border of the grid. Clearly if there is a realization of (m,n, h, v) then 

1. (vm, n, h, 1) and (m, hn, 1, v) must have a realization as well 

2. and every separated rectangle must admit a tiling with horizontal bars of length h and vertical bars of 
height u, even without any tomographic constraint. 

Left part of figure ^ shows an instance which satisfies the first but not the second condition, since each of 
the four separated rectangle has odd size. However the two conditions are not sufficient: The right hand 
instance satisfies conditions 1 and 2. But it has no solution, since the last column must be filled with vertical 
dominoes, the remaining cells of the first row must be tiled with horizontal dominoes, and for the remaining 
rectangle we end up with the left hand side instance. However the two conditions are not sufficient, as shows 
the right part. 
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